#include <stdio.h>
#include <stdlib.h>
#include <math.h>//引入sin x,cos x,e^x的库
//使用指向函数的指针变量来复用一个通用函数
int main()
{
 int i,j,k,n,m;
 float res;
 float (*fun)(float);//定义指向函数的指针变量
 float integral (float a,float b,float (*fun)(float),int n);//n为将积分区间（b-a)分成n等分，当n的值越大的时候结果越精确，数学上的定义是取无穷大
 float fsin(float);//sin x函数的声明
 float fcos(float);//cos x函数的声明
 float fexp(float);//e^x 函数的声明
 printf ("请输入积分的下限\n");
 scanf ("%d",&m);
 printf ("请输入积分的上限\n");
 scanf ("%d",&n);
 printf("请输入你要计算的函数的具体函数\n");
 printf("1.sin(x) 2.cos(x) 3.e^x\n");
 scanf ("%d",&i);
 switch(i)
 {
 case 1:
  fun=fsin;//函数地址（入口）交给指针变量，灵活性强
  break;
 
 case 2:
  fun=fcos;
  break;
 case 3:
  fun=fexp;
 }
 //res=(fun)(2.00);
 //printf("xxxx=%f\n",res);
 printf ("计算的结果为\n");
 res=integral(m,n,fun,20000);
 printf("res=%f",res);
 return 0;
}
float fsin(float x)
{
 // printf("fsinx=%f\n",x);  
 return sin(x);
}
float fcos(float x)
{
 //printf("fcosx=%f\n",x);
 return cos(x);
}
float fexp(float x)
{
 //getchar();
 //printf("fexp=%f\n",x);
 return exp(x);
}
float integral(float a,float b,float (*fun)(float),int n)
{
 //矩形法计算,定积分转换为连续求和的形式
 int i=0;
 float x=a,s=0;
 float h=(b-a)/n;
 
 for (i=1;i<=n;i++)
 {
  x=x+h;
  s=s+((*fun)(x)*h);
 
 }
 return s;
}